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ABSTRACT 

Aims. We model the present-day population of classical low-mass X-ray binaries (LMXBs) with neutron star accretors, which have 
hydrogen-rich donor stars. Their population is compared with that of hydrogen-deficient LMXBs, known as ultracompact X-ray 
binaries (UCXBs). We model the observable LMXB population and compare it to observations. We model the Galactic Bulge because 
it contains a well-observed population and it is the target of the Galactic Bulge Survey. 

Methods. We combine the binary population synthesis code SeBa with detailed LMXB evolutionary tracks to model the size and 
properties of the present-day LMXB population in the Galactic Bulge. Whether sources are persistent or transient, and what their 
instantaneous X-ray luminosities are, is predicted using the thermal-viscous disk instability model. 

Results. We find a population of ~2.1 x 10 3 LMXBs with neutron star accretors. Of these about 15 - 40 are expected to be persistent 
(depending on model assumptions), with luminosities higher than 10 35 erg s” 1 . About 7-20 transient sources are expected to be in 
outburst at any given time. Within a factor of two these numbers are consistent with the observed population of bright LMXBs in the 
Bulge. This gives credence to our prediction of the existence of a population of ~ 1.6 x 10 3 LMXBs with low donor masses that have 
gone through the period minimum, and have present-day mass transfer rates below 1(L H M s yr _1 . 

Conclusions. Even though the observed population of hydrogen-rich LMXBs in the Bulge is larger than the observed population of 
(hydrogen-deficient) UCXBs, the latter have a higher formation rate. While UCXBs may dominate the total LMXB population at the 
present time, the majority would be very faint or may have become detached and produced millisecond radio pulsars. In that case 
UCXBs would contribute significantly more to the formation of millisecond radio pulsars than hydrogen-rich LMXBs. 

Key words, binaries: close - Galaxy: bulge - X-rays: binaries 


1. Introduction 


Low-mass X-ray binaries (LMXBs) are binaries in which a star 
with a mass below ~ 1.5 AT. transfers mass via Roche-lobe over¬ 


flow to a neutron star or black hole companion (e.g., van den 
HeuveT] |T975[ [Joss & Rappaport||l979| l. The bifurcation period 
separates zero-age main-sequence LMXBs with expanding or¬ 
bits and long orbital periods, in which mass transfer is driven by 
donor evolution, from LMXBs with shrinking orbits, in which 
mass transfer is driven by angular momentum loss via magnetic 
braking and gravitational wave radiation ( Tutukov et al.|[T985) 
Pylyser & Savonije|1988|[l989| l. In the first group of ‘diverging’, 
long-period systems, the donor star will eventually become a de¬ 
tached white dwarf after its subgiant or giant progenitor has lost 
almost all of its envelope fWebbink et al.|1983| . In the ‘converg¬ 
ing’ systems, however, the donor loses most of its mass before 
it has developed a helium core. Convection turns the donor into 
a homogeneous mixture of hydrogen and helium, where the ra¬ 
tio between both elements depends on the evolutionary stage at 


the time hydrogen burning was extinguished (e.g., Tutukov et al. 
[1985) . 

Similar to cataclysmic variables, the orbital periods of the 
bulk of the converging LMXBs shrink until a period minimum 
of ~ 70 - 80 min (Paczynski 1981 1 Paczynski & Sienkiewicz 
|1981[ |Rappaport et al.||1982|>, where the donor reaches a max¬ 


imum average density, and later it becomes degenerate. In the 
case of sufficiently evolved cores, the period may continue to 
shrin k to ~40 min (|Nelson & Rappaport 20 03]) or even ~5 — 10 
min (|Tutukov et al.||1985[ |Podsiadlowski et aL]|2002[ |van der| 
|Sluys et al.|2005 ). Systems subsequently slowly expand towards 
longer orbital periods. Stable mass transfer continues to decrease 
the donor mass, and irradiation of the donor by the accretor and 
the accretion disk becomes important in driving mass transfer 
and generating a wind from the donor (Ruderman et al.|1989b|i. 


In both the converging and diverging systems, neutron star 
accretors can be recycled to spin periods of a few millisec¬ 
onds. If accretion stops, systems can become binary millisec¬ 
ond radio pulsars, as has been observed in the 4.75 hr binary 
FIRST J102347.67 +003841.2 ( jArchibald et al.l[2009] l. The ob¬ 
served evaporation process suggests the possibility that ‘black 
widow’ systems may turn into isolated millisecond radio pulsars 
( |Ruderman et al.|1989aHFruchter et al.|1988| l. 


In this study we estimate the number of classical (i.e., 
hydrogen-rich) LMXBs in the Galactic Bulge, and compare 
the results to observations of bright LMXBs, as well as with 
earlier population synthesis studies. Furthermore, we compare 
the results with the population of ultracompact X-ray binaries 
(UCXBs) in the Bulge modeled by |van Haaften et al7] ( |2013| l, 
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and to results by the Galactic Bulge Survey (Jonker et al. 2011[ 
[20l4l >. 


2. Method 


We use the binary population synthesis code SeBa (|Portegies 


Zwar t & Verbunt||19961 |Nelemans et al.|[2.Q01 1 |Toonen et al. 


2012[ > to simulate the primordial Galactic Bulge binary popula¬ 
tion, and select all binaries containing a neutron star and a main- 
sequence star, immediately after the supernova event. We use a 
combined common-envelope parameter o-oyl = 2 for massive 
stars (e.g., |Portegies Zwart & Yungelson 1998) 1. In Sect. [4] we 
also consider the effect of using a lower value. The initial mass 
function of the primary components is derived by combining the 
initial mass function for stellar systems ( Kroupaj2001 i with the 
binary fraction as a function of mass. The mass ratios of the sec¬ 
ondary and primary components follow a uniform distribution 
between 0 and 1, and the eccentricity e is drawn from a distri¬ 
bution proportional to e between 0 and 1. The semi-major axis 
a distribution is inversely proportional to a (Popova e t al.|1982| ), 
with a lower lim it determined by the stellar r adii, and an upper 
limit of 10 6 R a (iDuquennoy & Mayor 1991). The kick veloci¬ 
ties imparted on neutron stars during their formation are drawn 
from the distribution by PaczynskT] ( 1990| ) with a dispersion of 
270 km s. For more details on the initial binary parameters, 
we refer to |van Haaften et ak| ( |2013l >. 

For the subsequent evolution of the selected systems we use 
the low-mass X-ray binary tracks described in |van der S luys 
et al. | ( |2005 1, a selection of which is shown in Fig. [T] These 
tracks describe the evolution of binaries starting with a neutron 
star accretor and a zero-age main-sequence donor. The tracks 
cover a grid of initial donor masses and orbital periods. The 
donor masses range from 0.7 to 1.5 M 0 with steps of 0.1 M 0 . 
The initial orbital periods range from 0.5 to 2.75 days, with steps 
of 0.25 days, but with 0.05 day steps near the bifurcation period. 
The spikes in Fig. [T| have a numerical origin; they are caused 
by the finite number of grid points in the stellar-structure mod¬ 
els. When a convective region shrinks or expands, a grid point 
near the boundary will start or stop being convective. Hence, 
the changes in convective regions are not smooth, but undergo 
discrete steps. Such changes cause steps in the radius of the star, 
which influences its Roche-lobe-filling factor, to which the mass- 
transfer rate is rather sensitive, causing the spikes. Each system 
produced by SeBa is matched to the track that most closely ap¬ 
proaches its donor mass and orbital period. We correct the orbital 
period of the SeBa system because its donor is slightly evolved 
at the time of the supernova event, whereas the tracks start with 
zero-age main-sequence donors. 

In Fig. □ the orbital periods have shrunk before the systems 
start mass transfer. The two tracks on the right-hand side show 
LMXBs that diverge, and ultimately become detached. The three 
tracks on the left-hand side represent converging systems, which 
go through a period minimum. Among the converging systems, 
the lowest period minimums are reached by systems with the 
longest initial period, as long it lies below the bifurcation period. 

Initially the orbital evolution of the binary is predominantly 
driven by magnetic braking of the main-sequence star. We calcu¬ 
late the angular momentum loss rate using the formula by |Rap-| 
paport et al.]( 1983| >. A metallicity Z = 0.02 is used in SeBa. The 
LMXB tracks use Z = 0.01, but their evolution is very similar 
for Z = 0.02 (van der Sluys et al.|2005j). 

P0l3) ,' 


Following van Haaften et al. (|2013|i, we approximate the star 


formation history of the Galactic Bulge by a Gaussian distribu¬ 
tion with a mean p — -10 Gyr (the present time is defined as 0) 
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Fig. 1 . Sample of tracks showing the mass transfer rate against orbital 
period for LMXBs with initial donor masses of 1 M Q . The complete 
grid covers a range of initial orbital periods between 0.5 and 2.75 days 
(12 - 66 hr, indicated by the horizontal bar) {van der Sluys et al.|2005|l. 


and a standard deviation <x of either 0.5 or 2.5 Gyr, in order to 
model both burst-like and extended star formation. For the total 
star forming mass we take 1 x 10 10 M Q 
|Wyse[2009 j). 

The X-ray luminosity Lx of a system is calculated from its 
mass transfer rate using 


(Clarkson & Rich 2009 


Lx ~ 7bol“ 


GM a M a 


R a 


(1) 


where r/b„i ~ 0.55 is the bolometric correction (following |Bel-| 
|czynski et al.|2008] >, M a is the accretor mass, R. d is the accretor 
radius, and G the gravitational constant. We use the thermal- 
viscous disk instability model (see, e.g., the review by |Lasota] 
|2001j ) to decide which of our modeled sources are transient and 
which are persistent at the present time. Sources do not have 
a stable disk if their mass transfer rates are below the critical 
value M cr i t for irradiated accretion disks derived by 
( 1999| , in the form given by in’t Zand et al. ( 2007| 


Dubus et al. 


M C1 


5.3 x KT 11 / 


M, 


0.3 


M., \ IP or b \ 


1.4 


V hr / 


M g yr 


( 2 ) 


where P or b is the orbital period and / is a scale factor depending 
on the disk composition; / « 1 for solar-composition disks and 
f ~ 6 for helium disks. Systems with unstable accretion disks 
are expected to be transient sources, which show outbursts sepa¬ 
rated by (much longer) quiescent phases. The luminosity during 
outbursts is given by Eq. (Jj]l with accretion rate M a = M cr ;t- 
The luminosity is assumed to be zero during quiescence. There¬ 
fore, the duty cycle is equal to the ratio of the time-averaged 
mass transfer rate and the mass transfer rate during outburst. The 
time-averaged mass transfer rate is the rate at which the donor 
loses mass according to the evolutionary tracks. This rate is vir¬ 
tually constant on timescales much shorter than the evolution¬ 
ary timescales of the donor and the binary system. If sources 
have time-averaged mass transfer rates higher than the critical 
value, we assume they are persistent X-ray sources with X-ray 
luminosities given by Eq. (|T]» where M a is taken to be the time- 
averaged mass transfer rate. 

We do not model LMXBs with black hole accretors because 
it is not at all clear how they are formed (Portegies Zwart et al. 
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Fig. 2. Model distribution of present-day mass transfer rate versus or¬ 
bital period for hydrogen-rich LMXBs in the Galactic Bulge. White 
squares correspond to 0.4 or fewer LMXBs. The highest density in the 
figure is 10 27 ~ 500 per pixel. The black line shows the critical mass 
transfer rate for thermal-viscous disk instability, given by Eq. 0 with 
M a = 1.4 M 0 and / = 1 (solar composition). The star formation history 
width cr = 0.5 Gyr. 


1997 

Kalogera |1999( Podsiadlowski et al. 2003 

Justham et al.| 

2006 

Yungelson et al. 2006 ), and in particular 

tow the binary 


system remains bound during the formation of the black hole, 
given the low mass of the companion. Additionally, all observed 
(certain) black hole LMXBs are transient and therefore we can¬ 
not compare model results to persistent systems. 


3. Results 


3. 1. The total population of LMXBs 


The total number of LMXBs with neutron star accretors formed 
over the history of the Galactic Bulge in our study is 7.6x 10 3 . Of 
these, (2.0 - 2.1) x 10 3 (for star formation history width cr = 2.5 
and 0.5 Gyr, respectively) are still X-ray binaries at the present 
time - all other systems have stopped transferring mass and have 
turned into long-period detached white dwarf-neutron star bina¬ 
ries. 


Figures[2]and[3]show the present-day population of ~2 x 10 3 
LMXBs in terms of their orbital periods and time-averaged mass 
transfer rates, for the two star formation history widths (Sect. |2|». 
There are (1.5 - 1.6) x 10 3 systems (for cr = 2.5 and 0.5 Gyr, 
respectively), three quarters of all LMXBs, with mass transfer 
rates below 10 1 M G yr 1 and orbital periods shorter than 3 hr. 
These LMXBs have passed through the orbital period minimum 
and have increasing periods at the present time. The donor stars 
in these systems are (semi-)degenerate and have masses below 
~0.05 M q . The donor stars in systems that reach periods shorter 
than ~ 60 min have no hydrogen in their cores and less than 10% 
hydrogen on their surface at the period minimum (van der Sluys 
let al.|20051 >. 


Almost all LMXBs with orbital periods longer than ~ 3 hr 
have mass transfer rates above ~ 10 _l1 M Q yr -1 , and about three 
quarters of these are diverging systems. The solid lines in Figs. [2] 
and [3] separate the persistent and transient sources according to 
the Disk Instability Model (Sect. [2]). In Sect. 3.2 we will discuss 
this in more detail. 


Fig. 3. Same as Fig. [5] except cr = 2.5 Gyr. The highest density in the 
figure is 10 24 « 250 per pixel. 



Fig. 4. Orbital period distributions of the present-day populations of 
hydrogen-rich LMXBs with neutron stars accretors (solid line), and 
LMXBs that have evolved into detached white dwarf-neutron star sys¬ 
tems (dashed line) in the Galactic Bulge. The star formation history 
width cr = 0.5 Gyr. 


The orbital period distributions of all present-day neutron 
star systems, both semi-detached and detached, are shown in 
Figs. [4] and [5] again for two star formation histories. The semi¬ 
detached systems (solid lines) are the same populations as shown 
in Figs.[2]and[3j combining all mass transferrates. Of all LMXBs 
that have formed in the history of the Bulge, about 78% are 
diverging systems. The majority (~ 80%) of present-day semi¬ 
detached systems (solid lines) have descended from converg¬ 
ing LMXBs. All modeled present-day detached systems (dashed 
lines) have descended from diverging LMXBs, of which the 
donor envelope has been lost and the core has been left as a 
helium or carbon-oxygen white dwarf. The orbital periods of 
the white dwarf-neutron star systems are ~ 6 hr or longer. The 
shortest-period detached systems may have turned into UCXBs, 
if they have had enough time to shrink their orbit, and if they re¬ 


mained stable during the onset of mass transfer (Yungelson et al. 
|2002||van Haaften et al.|2012b| l. 
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Fig. 6. Time-averaged mass transfer rate distributions of the present-day 
population of hydrogen-rich LMXBs with neutron star accretors in the 
Galactic Bulge. The star formation history width is cr = 0.5 Gyr (solid 
line) and cr = 2.5 Gyr (dotted line). 

The orbital periods of the LMXBs vary between about 20 
min and 10 d, with a peak around 2.5 hr. Depending on donor 
mass, the bifurcation period is in the range 0.8 - 1.5 d. The long- 
period end of the LMXB distribution represents diverging sys¬ 
tems. The donor stars in present-day LMXBs in the Bulge origi¬ 
nally had masses lower than ~ 1 M 0 . 

Figure [6] shows the time-averaged mass transfer rates of the 
neutron star LMXB population. About 330 systems have time- 
averaged mass transferrates exceeding 10 111 M 0 yr 1 , approx¬ 
imately 1% of the Eddington limit. The post-period minimum 
systems have lower average mass transfer rates than the popula¬ 
tion of diverging, long-period systems. 

3.2. The observable population of LMXBs 

As discussed in Sect. [2j only the brightest and shortest-period 
sources in the modeled population are persistent, while the re¬ 
mainder of the population experience outbursts, with a duty cy¬ 
cle that decreases towards longer orbital period and lower av¬ 
erage mass transfer rate. At the present time, the number of 



Fig. 7. Orbital period distributions of the present-day population of per¬ 
sistent hydrogen-rich LMXBs with neutron star accretors in the Galactic 
Bulge (black lines), as well as a snapshot in time of the transient popula¬ 
tion in outburst (gray lines). The star formation history width is cr = 0.5 
Gyr (solid lines) and cr = 2.5 Gyr (dashed lines). Unlike the previous 
histograms, here the vertical scale is linear. 

LMXBs in the persistent regime is 37 for cr = 0.5 Gyr and 42 
for cr = 2.5 Gyr, whereas the number of transient LMXBs in out¬ 
burst is 15 for cr = 0.5 Gyr and 20 for cr = 2.5 Gyr. Figures[7 and 
[8] show that almost all persistent systems have orbital periods in 
the range 0.3 - 10 hr and X-ray luminosities of 1 () '"’ 37,3 erg s 1 . 
The transient sources in outburst are brighter on average than the 
persistent sources. The orbital periods of transients in outburst 
lie in the range 0.5 - 30 hr, with a peak near 2 hr, and their lumi¬ 
nosities lie between 10 36 erg s 1 and the Eddington luminosity. 
The present-day observable population, both the persistent and 
transient sources, are predominantly of the converging type, es¬ 
pecially for a narrow star formation history (cr = 0.5 Gyr). The 
transient sources have mostly evolved beyond the period mini¬ 
mum. These correspond to the darkest squares in Figs. [2] and [3] 
The persistent sources are typically near, or on their way to, the 
period minimum. 

It is interesting to note that the period distribution of the 
bright, observable population (Fig. [7j does not show a global 
maximum near the period minimum. Instead, the distribution is 
relatively uniform between 20 min and 10 hr. This is similar to 
the situation for Cataclysmic Variables, where a peak near the 
period minimum appears only after including very faint sources 
( |Gansicke et al.||2009| l. In the LMXB case, our models predict 
an accumulation of semi-detached systems beyond the period 
minimum near 2 hr (solid lines in Figs.[4]an d0, where angular 
momentum loss via gravitational wave radiation becomes very 
weak, and evolutionary timescales long. However, the vast ma¬ 
jority of these systems is transient as a result of their low mass 
transfer rates. Because of their low duty cycles, only a small frac¬ 
tion is in outburst at a given time. On the other hand, systems 
with periods near 3 - 10 hr are either persistent, or transient with 
relatively high duty cycles. Therefore there is no peak in the ob¬ 
servable population at short orbital periods. 


4. Discussion 

We predict the existence of ~ 40 persistent LMXBs and, at any 
given time, ~ 20 transient LMXBs in outburst in the Galactic 
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Fig. 8. X-ray luminosity distributions of the present-day population 
of persistent hydrogen-rich LMXBs with neutron star accretors in the 
Galactic Bulge (black lines), where mass transfer rates have been con¬ 
verted to X-ray luminosities using Eq. 0. as well as a snapshot in time 
of the transient population in outburst (gray lines). The star formation 
history width is a = 0.5 Gyr (solid lines) and <x = 2.5 Gyr (dashed 
lines). 


Bulge. These numbers are not very sensitive to the shape of the 
star formation history function. 

However, the numbers of predicted persistent LMXBs 
(Sect. [3T2] ) are uncertain for several reasons. There are uncertain¬ 
ties in the population synthesis and LMXB tracks, but also in the 
disk irradiation (Dubus et al. 1999). The latter is important be¬ 
cause the critical lines in Figs. [2] and [3] coincide with relatively 
narrow bands of LMXBs. For instance, increasing the critical 
mass transfer rate by a factor of 2 would reduce the number of 
persistent sources by a factor of 2 to 4. A different critical mass 
transfer rate would also affect the duty cycle of transient sources. 


If we use a lower value for the combined common-envelope 
parameter ac+A = 0.2 in SeBa, rather than the value of 2 
(Sect. |2j, as well as a Maxwellian velocity distribution of the 
neutron star kicks with a dispersion of 450 km s , the total num¬ 
ber of LMXBs formed is smaller by a factor of ~ 2.7. The rel¬ 
ative numbers of persistent and transient systems do not change 
significantly. 

Because the masses of neutron stars at birth are very un¬ 
certain, we adopted an initial neutron star mass of 1.4 M 0 in 
our simulations. However,[Kuranov et al.|( |2014j ) recently found 
that a uniform distribution of initial neutron star masses between 
1.4- 1.9 M 0 gives the best match with observations of persis¬ 
tent LMXBs in the Galactic Bulge. More massive neutron stars 
would affect our results in several ways. A more massive super¬ 
nova remnant implies a higher probability for the binary system 
to survive the supernova event during which the neutron star is 
formed. For a given Roche-lobe filling donor star, a higher ac- 
cretor mass will cause a higher time-averaged mass transfer rate 
if mass transfer is driven by angular momentum loss via grav - 
itational wave radiation (M oc M[ 


2/3 


van Haaften et 


al. 2012b). 


Additionally, a more massive accretor and a higher mass transfer 
rate both directly lead to a higher average luminosity (see Eq.[T| 
where the neutron star mass-radius relation should also be taken 
into account). The critical mass transfer rate (Eq.[2| also depends 
on accretor mass, although rather weakly (oc Thus, the ef¬ 

fect of a more massive accretor is as follows. Systems that are 


already persistent will become brighter. Transient systems close 
to the critical mass transfer rate may become persistent, however, 
persistent systems in which angular momentum loss is caused by 
magnetic braking may become transient. Transient systems well 
below the critical line will remain transient, but have a higher 
duty cycle. A larger spread in neutron star masses will lead to a 
larger spread in observed X-ray luminosities as well. 

Several studies ( |Pfahl et al.|2002||Podsiadlowski et al.|2004| ) 
have suggested the existence of a dichotomy in neutron star 
kicks, where the kick velocity (and neutron star mass) depends 
on the type of supernova the neutron star is formed in, which in 
turn depends on the degree to which the primary was influenced 
by the companion during its evolution. In our standard simu¬ 
lation, about 43% of the neutron stars that end up in LMXBs 
originated from 8-11 M 0 primaries, which is a necessary con¬ 
dition for the formation of an electron capture supernova ( |Podsi-| 
|adlowski et al.|2 004 ). Neutron stars forming in these supernovae 
are expected to have very low kick velocities (< 50 km s 1 ). 
Out of all neutron stars in LMXBs in our simulation, about 26% 
went through a helium star/helium giant phase where the initial 
helium star mass was in the range of 2.6 - 2.95 M 0 , which ac¬ 
cording to calculations by Tauris et al. ( |2015| ) can lead to an elec¬ 
tron capture supernova, depending on the initial orbital period. 
The Paczynski] (1990) kick velocity distribution, which we use 
in our standard model, includes a sufficient fraction of low kick 
velocities to account for the neutron stars formed during elec¬ 
tron capture supernovae. Another reason we use the Paczynski 
kick distribution is to be able to directly compare the results to 
our earlier work on UCXBs ( |van Haaften et al.|2013| ), where we 
used the same distribution. 

Several population synthesis studies of LMXBs have been 
performed in the past. Fragos et ah] ( j2008) > predicted ~ 800 
LMXBs with neutron star accretors with Lx > 10 36 erg s in 
an elliptical galaxy (with an age similar to the Bulge), most of 
which are persistent. The stellar mass of the modeled elliptical 
galaxy is nine times larger than the stellar mass we use for the 
Bulge, therefore their estimate is equivalent to ~ 90 LMXBs in 
the Bulge. We find a total of 23 - 31 persistent LMXBs with lu¬ 
minosities above this value in the Bulge, and at any given time 
another 15 - 20 in outburst (for cr = 0.5 and 2.5 Gyr, respec¬ 
tively), see Fig. [8] Hence, our number of sources that is brighter 
than 10 36 erg s ' is a factor of ~ 2 - 3 smaller, per unit stellar 
mass. Fragos et al.|( |2008| ) found most LMXBs near a period of 
~ 1 hr, with a lower peak near 30 hr. This resembles the dis¬ 
tribution in our Figs. [4] and [5] (solid lines), which show peaks 
near 2 and 15 hr. |Kalogera & Webbink| ( |1998) > and jKalogera] 
( |1998| ) found a formation rate of LMXBs with neutron star ac¬ 
cretors of 10 6 — 10 s yr 1 in the Galaxy, depending on common- 
envelope parameters and supernova kick velocities. This implies 
a total number of LMXBs formed over the history of the Disk of 
roughly 10 4 ~ 5 , which compares reasonably well to our result of 
~ 10 4 , taking into account that the Bulge has a stellar mass that is 
4-6 times lower than the stellar mass of the Disk (Klypin et al. 


2002|). Kiel & Hurley ( 2006|), in their favored model, found about 


1700 LMXBs with neutron star accretors in the Galaxy. |Kiel| 
|& Taam| ^2013) ) modeled the population of LMXBs and black 
widow systems. Their various models produce 5 - 200 black 
widow systems, and about one hundred times as many LMXBs 
in the Galaxy. These numbers are consistent with our results 
as well, again after correcting for the mass difference between 
Bulge and Galaxy. Kuranov e t ak| (2014) modeled an old stellar 
population of 500 000 systems. This corresponds to a total stellar 
mass of approximately 5 x 10 7 - 10 8 M 0 , based on the stellar ini¬ 
tial mass function. Therefore the numbers these authors predict. 

Article number, page 5 of[8] 



































































A&A proofs: manuscript no. vanhaaften_etal_2015_lmxb 


about 0.1-1 LMXBs in their population across their models, 
should be multiplied by 100 - 200 in order to compare them 
to our numbers for the Galactic Bulge. The number of LMXBs 
they found adjusted this way is a factor of 0.3 - 2 times our num¬ 
ber, depending on the models used. Their LMXBs are generally 
somewhat brighter than the systems we modeled, in particular 
they find more systems with luminosities over 10 37 erg s 1 . We 
note that agreement between different models within a factor of 
a few is still acceptable, given for instance the uncertainty of a 
factor of ~ 3 based on the common-envelope parameters and kick 
velocity distribution. 

The model results must be compared to the observed pop¬ 
ulation of LMXBs, in particular bright systems. In general this 
is difficult because of the variability of LMXBs. Voss & Ajello 
( 2010 1 found 18 LMXBs in hard X-rays (15-55 keV) in data 
from the Swift Burst Alert Telescope, within 10° (1.5 kpc) of 
the Galactic Center, that are not known to be located outside 
the Bulge (e.g., in a globular cluster). Owing to confusion with 
background sources near the Galactic Center, there may be more 
sources. On the other hand, as a result of poor distance con¬ 
straints, some of the observed sources may be located in front 
of, or behind, the Bulge. The faintest of the 18 LMXBs has a lu¬ 
minosity of ~ 2.4x 10 35 erg s 1 in the hard X-ray band, which im¬ 
plies a total X-ray luminosity of ~ 7x 10 35 erg s _1 after a spectral 
correction based on Revnivtsev e t allj ( {2008 ). The number of 18 
can be compared to our prediction in Fig. [8 which shows about 
45 sources above this luminosity threshold. The LMXB catalog 
by Liu et al. ( 2007) contains approximately 40-50 LMXBs (per¬ 
sistent and transient sources combined) that are located in the 
direction of the Bulge based on their celestial coordinates. Fur¬ 
thermore these sources are not known to be located in globular 
clusters, neither are they candidates of having black hole accre- 
tors. About one-third of these sources have a distance estimate 
that is consistent with the distance to, and size of, the Bulge. The 
others are located sufficiently close to the Galactic Center (e.g., 
60% within two degrees) that it is probable that the vast major¬ 
ity of them is located in the Bulge, given the high stellar density 
near the Galactic Center. The Liu et al.] (2007 1 catalog includes 
sources found by pointed observations, therefore the sample is 
inhomogeneous. These numbers lie within a factor of two of our 
prediction of bright X-ray sources in the Bulge. A lower value 
for ffcE'k as discussed above, would give an even better match. 
We consider the agreement of the number of predicted systems 
with the observed ones good. We compare the orbital periods 
and luminosities of our modeled systems to persistent Galactic 
LMXBs in Revnivtsev et alJ( 201l) . The X-ray luminosities of 
those systems in the 2-10 keV range are ~ io 36 - 5-38 erg s _1 at or¬ 
bital periods of ~ 2-5 hr. This corresponds well with the proper¬ 
ties of our persistent LMXBs, although we also find fainter sys¬ 
tems with luminosities of io 35 5-36 5 erg s 1 . The Galactic Bulge 
Survey, which mostly aims to detect quiescent LMXBs ( |Jonker| 
et al .|2011 |2014[ >, has observed 12 square degrees in the Bulge. 
Extrapolating from our prediction of ~ 60 bright hydrogen-rich 

50 UCXBs 


LMXBs in the Bulge as well as 


van Haaften et al. 


2013)), and assuming they are evenly distributed over the pro 


jected area of the Bulge, this implies up to six bright sources in 
the Galactic Bulge Survey field. 


4.1. Comparison with the population of ultracompact X-ray 
binaries 

Ultracompact X-ray binaries (UCXBs) are hydrogen-deficient 
LMXBs, with compact donors and observed orbital periods 
shorter than 1 hr (e.g., [Nelemans & Jonker pOlO) |van Haaften | 
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Fig. 9. Model distribution of present-day mass transfer rate versus or¬ 
bital period for all UCXBs (the vast majority of systems to the left of the 
dotted line) and hydrogen-rich LMXBs (virtually all systems to the right 
of the dotted line) in the Galactic Bulge, assuming the UCXB evolution 
is driven only by gravitational wave radiation, and that all systems re¬ 
main semi-detached. White squares correspond to one or fewer systems. 
The highest density in the figure is 10 5 (outside the scale). The solid line 
shows the critical mass transfer rate for thermal-viscous disk instability, 
given by Eq. 0 with M. d = 1.4 M 0 and / = 6 (helium composition), 
and likewise, the dashed line for / = 1 (solar composition). The star 
formation history width <x = 0.5 Gyr. The hydrogen-rich LMXB popu¬ 
lation in this figure is the same as in Fig. [2] but a different density scale 
is used. The UCXB data in this figure are taken fromlvan Haaften et al.l 
(20T3l l. 


|et al.|[2012b| |Heinke et al.||20f3] l. As mentioned in Sect. [T] 
UCXBs can also form via an evolved main-sequence donor, but 
more than 99% of all UCXBs is expected to form after a sec¬ 
ond common envelope phase, when eventually the core of the 
secondary fills its Roche lobe (van Haaften et al. 2013) . 

Figure 9] shows the present-day UCXB population in the 
Bulge (to the left of the dotted line) assuming their evolution is 
not affected by an irradiation-induced wind from the donor (van 
[Haaften et al.||2013]>, as well as the population of hydrogen-rich 
LMXBs from Fig. [2](to the right of the dotted line) for compar¬ 
ison. The oldest UCXBs have the lowest mass transfer rates and 
the longest orbital periods. Without a donor wind, the evolution 
of an UCXB is driven by angular momentum loss only via grav¬ 
itational wave radiation, and the system remains semi-detached. 
In the case of a wind from the donor, the mass transfer rate is 
higher, and longer orbital periods can be reached within the age 
of the Universe ( |van Haaften et al.|2012a| l. Donor evaporation is 
expected to have an important impact on the orbital period dis¬ 
tribution, which will shift towards longer periods and reduce the 
number of low mass transfer rate systems. 

Haaften et al. (2013|) predicted a population of (0.2 - 


1.9) X 1(F UCXBs in the Bulge, depending on model assump¬ 
tions, of which a large fraction may have turned into (binary) 
millisecond pulsars as a result of recycling and accretion turnoff. 
When we compare the population of UCXBs in Fig. ^(show¬ 
ing 1.9 x 10 5 UCXBs) to the hydrogen-rich LMXBs (2.1 x 10 3 
systems), we see that UCXBs are about one hundred times more 
common than hydrogen-rich LMXBs, unless the vast majority 
of UCXBs have become detached and turned into millisecond 
radio pulsars (for the lower value of a ce 4 = 0.2 the ratio of 
~ 100 decreases to ~ 30). In the case UCXBs typically become 
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detached, hydrogen-rich LMXBs, such as the evaporating black 
widow systems, cannot be the main progenitor class of (iso¬ 
lated) millisecond radio pulsars. In fact, if most UCXBs become 
isolated millisecond radio pulsars, their absolute number corre¬ 
sponds reasonably well with the estimated number of isolated 
millisecond radio pulsars (van Haaften et al. |2013| l. Therefore, 
apart from the relative importance of UCXBs, the absolute num¬ 
bers also suggest that the classical LMXBs cannot explain the 
size of the millisecond pulsar population. The higher birth rate 
of UCXBs compared to hydrogen-rich LMXBs is probably a re¬ 
sult of the higher secondary masses in zero-age main-sequence 
UCXB progenitor binaries, which make it easier to avoid disrup¬ 
tion of the system during the supernova explosion of the primary, 
and to survive the common envelope stage. 

If a significant fraction of UCXBs do not become millisec¬ 
ond radio pulsars, their number, (0.2 - 1.9) x 10 5 , is much higher 
than the number of hydrogen-rich LMXBs, (0.8 - 2) x 10 3 (for 
the same variation in model parameters). However, the latter are 
expected to be slightly more common in the bright population 


(> 10 35 erg s '), with 20 - 60 versus 5-50 for UCXBs (van 


|Haaften et al.|20 13). Thus, a much larger fraction of hydrogen- 
rich LMXBs is bright, compared to UCXBs. This has several 
reasons. The donor composition in UCXBs can be helium, car¬ 
bon and oxygen, or a mixture of these elements. LMXBs that 
start mass transfer from a hydrogen-rich donor have donors com¬ 
posed of hydrogen and/or helium. At low donor mass, the donors 
composed of lighter elements are larger than donors of heav¬ 
ier elements (of the same mass), due to Coulomb interactions 
( Zapolsky & Salpeter|1969 1. Therefore, the donors composed of 
heavier elements reach the same average donor density (hence 
also the same orbital period, by the period-density relation) at a 
lower donor mass. Because gravitational wave losses are smaller 
at lower donor mass, UCXBs have a lower time-averaged mass 
transfer rate at the same orbital period (and average density) 
than LMXBs with hydrogen in their donors (which can be seen 
in Fig. [9}. Moreover, the critical mass transfer rates are higher 
for UCXBs (solid and dashed line in Fig. [9] respectively). Com¬ 
bined, this leads to significantly lower duty cycles for UCXBs - 
this explains why the number of UCXBs in outburst is not larger 
than that of hydrogen-rich LMXBs in outburst, even though the 
UCXBs are much more common (unless they have become mil¬ 
lisecond pulsars). Persistent sources have a higher critical mass 
transfer rate that results in a shorter life as a persistent source. 
The result is that a larger fraction of the population of hydrogen- 
rich LMXBs is persistent. 

In Fig. [TO] we combine the luminosity distribution of 
hydrogen-rich LMXBs (Fig.|8]l with the luminosity distribution 
of UCXBs by van Haaften et ak| ( |2013) >, both the persistent 
sources and the transient sources in outburst, for our standard 
model (cuceA = 2). In this model, the contributions of persistent 
and transient systems are approximately equal, but in the case of 
a common envelope that is ten times less efficient, the number 
of UCXBs decreases by a factor of ~ 8, whereas the number of 
hydrogen-rich LMXBs decreases by a factor of ~2.1. Therefore 
the latter are more common among bright sources by a factor of 
a few with this change in assumptions. 

The bright end of our luminosity distribution (Fig. [IQ]) can be 
compared with extragalactic LMXB observations parametrized 
by |Gilfanov| ( |2004| , shown as the dotted curve in Fig. [TO] The 
absolute numbers correspond quite well with our two models 
(solid and dashed bars). Counting only sources with luminosi¬ 
ties between 10 35 - 10 38 5 erg s ' 1 , our cr = 0.5 Gyr model has 90 
sources, and our cr = 2.5 Gyr model 110. Gilfanov (20041 has 
120 sources in this luminosity range. Our alternative ach^ = 0.2 



Fig. 10. X-ray luminosity distributions of the present-day population of 
LMXBs, including UCXBs, with neutron star accretors in the Galactic 
Bulge. The star formation history width is cr = 0.5 Gyr (solid line) and 
cr = 2.5 Gyr (dashed line). These distributions include hydrogen-rich 
as well as hydrogen-poor LMXBs, and both persistent sources and tran¬ 
sient sources in outburst. The dotted line shows the X-ray luminosity 
function parametrized by |Gilfanov[ j2004] sample ‘AH' in Table 3) us¬ 
ing the best-fitting normalization of K\ = 440.4, and further adjusted 
for the bin width of our histograms to allow for direct comparison. 


model produces about five times fewer sources than Gilfanov 


( |2Q04) >. The bright-end cut-off is located at a somewhat lower 
luminosity in our models, and is more gradual. Various uncer¬ 
tainties in the models could cause this. For instance, adopting a 
higher outburst luminosity for transient sources would make the 
distributions more similar. The decrease at low luminosities in 
our distribution is a result of not modeling quiescent stages, and 
in the case of transient sources, of assuming an abrupt transition 
from quiescence to outburst and back. 

Figure [TO] also shows that our model including recent star 
formation (cr = 2.5 Gyr, histogram with dashed line) shows a 
stronger break near 10 37 erg s than the model without recent 
star formation (solid line). This suggests that particulars of the 
star formation history could be a factor in explaining why spiral 
galaxies like the Milky Way (Grimm et al. 2002ji and M 31 (Kong 


|et al.|2003]|Trudolyubov & Pr iedhorsky 20041 show a break near 
this luminosity, whereas this break is probably absent in the old 
elliptical galaxies NGC 3379 and NGC 4278 ( |Kim et al.|2006] l. 


5. Conclusions 

We have simulated the evolution of binaries in the Galactic 
Bulge in order to predict the present-day population of low-mass 
X-ray binaries and their possible descendants, such as detached 
white dwarf-neutron star systems and millisecond radio pulsars. 
We predict about 2.0 x 10 3 LMXBs with neutron star accretors 
in the Galactic Bulge. Based on the predicted number of per¬ 
sistent sources and transient sources in outburst, this number 
corresponds well with observations. Because our prediction is 
supported by observations combined with the timescale progres¬ 
sion along the evolutionary tracks for LMXBs (which are more 
accurate than population synthesis), it is not very sensitive to un¬ 
certainties in population synthesis models. About three-quarters 
of the LMXBs are predicted to be systems below the bifurcation 
period that have passed through the period minimum, and have 
very low mass transfer rates and orbital periods shorter than 3 
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hr, some of which may have become millisecond radio pulsars. 
Most LMXBs that have formed over the history of the Bulge are 
no longer transferring mass at the present time, and we predict 
about 5.5x 10 3 binary millisecond radio pulsars with white dwarf 
companions in the Bulge with a LMXB origin (i.e., excluding the 
contribution of intermediate-mass X-ray binaries). 

Ultracompact X-ray binaries have a significantly higher for¬ 
mation rate than hydrogen-rich LMXBs and are the most com¬ 
mon subclass of low-mass X-ray binaries, unless the majority 
have become detached and turned into millisecond radio pulsars, 
in which case they dominate the millisecond radio pulsar forma¬ 
tion from all LMXBs. Even so, hydrogen-rich LMXBs are more 
easily observed because a larger fraction of them is persistent 
and the transient sources have longer duty cycles. 
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